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Ultracompact stars with multiple necks 
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We discuss ultracompact stellar objects which have multiple necks in their optical geometry. 
There are in fact physically reasonable equations of state for which the number of necks can 
be arbitrarily large. The proofs of these statements rely on a recent regularized formulation 
of the field equations for static spherically symmetric models due to Nilsson and Uggla. We 
discuss in particular the equation of state p — p — p s which plays a central role in this context. 



The purpose of this communication is to dis- 
cuss the existence of ultracompact relativistic stellar 
models which have multiple necks in their optical 
geometries. A few examples of such models were 
given in JJJ. Abramowicz and coworkers discussed 
the optical geometry of uniform density stars in 0] . 
Such stars have a single neck and an accompanying 
bulge. Necks (and bulges) correspond to unstable 
(stable) circular null orbits. We say that a star is 
ultracompact if the spacetime admits one or more 
spatially closed null orbits. This definition of ultra- 
compactness was proposed in S . As is well-known, 
a Schwarzschild black hole has exactly one closed 
null orbit at r — 3M. A uniform density star is 
ultracompact if its radius satisfies R < 3M. How- 
ever, for more realistic equations of state, there is 
no such simple relation between radius and ultra- 
compactness J3J. In the past, ultracompactness was 
usually defined by the condition R < 2>M. However, 
as was shown in |3[] that definition does not capture 
what may be regarded as the essence of ultracom- 
pact bodies, namely their ability to serve as a grav- 
itational trap for relativistic matter such as neutri- 
nos, electromagnetic waves or gravitational waves. 
In particular, the neck (or necks) of a compact star 
can be located entirely within the stellar interior. 

The optical geometry is a very useful tool to vi- 
sualize the structure of the null geodesies of static 
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spherically symmetric (SSS) models. Models with 
R < 3M have a neck in the exterior optical ge- 
ometry which is associated with the Schwarzschild 
circular null orbit at R = ZM . In the interior of 
the star there are also other spatially closed null 
orbits. This was believed to be the typical situ- 
ation for ultracompact stars. However, as will be 
discussed below, the optical geometries for a large 
class of more realistic equations of state have multi- 
ple necks. The maximum number of necks, n max , 
for a given equation of state is sensitive to the 
nunrber 7^, = 1 + linip^oo p/ p (assuming that it 
exists). If the high density equation of state has 
the Zel'dovich stiff matter form, p = p, so that 
7oo = 2, then the number of necks becomes arbi- 
trarily large as the central density, p Cl tends to in- 
ifinity. On the other hand consider the equation of 
state p = (7 — l)(p — p s ) where 7(= 700) and p s 
are constants. In the limit 7 — > 00 it goes over into 
the uniform density equation of state for which we 
know that n max = 1. Thus n max varies in the range 
[1, 00) as 7 varies from infinity to 2. At the other end 
n max = for 7 below a certain value 70 > 4/3 Q. 
These results and others mentioned in this commu- 
nication will be proved and discussed in more detail 
in|. 

The above statements can be inferred from a pow- 
erful formulation of the SSS field equations due to 
Nilsson and Uggla B, [fl. The idea is to replace the 
Tolman- Oppenheimer- Volkoff (TOV) equations by a 
maximally reduced dynamical system which has a 



state space which is both regular and compact. In 
fact Nilsson and Uggla derives three such systems, 
one adapted to linear equations of state H , and two 
complementary for polytropic equations of state M . 
We will use the Y>-K-y system given in JJJ, to be 
referred to as the loaf of bread system and the loaf 
of bread state space due to the shape of the state 
space. A price one has to pay is that the system is 
necessarily 3-dimensional and so has one extra di- 
mension compared to the TOV system. As a result, 
the loaf of bread system incorporates models which 
are irregular as well as regular at the center. In fact 
the inclusion of the irregular models is an essential 
feature in order to understand the structure of the 
solution space of the regular models. The under- 
lying reason for this will be touched upon below. 
For simplicity, in this paper we focus attention on 
the stiff fluid equation of state p = p — p s where 
p s ^ and the subindex "s" stands for evaluation 
at the surface of the star where p — 0. This is the 
most general equation of state for which the speed 
of sound (defined as ^Jdp/dp) equals the speed of 
light. In particular, there is no manifest violation of 
causality in contrast to the case of uniform density 
matter. 

We now proceed to show how the loaf of bread 
dynamical system and state space for SSS stars can 
be used to analyse compactness characteristics such 
as the occurrence of spatially closed null orbits. We 
begin by introducing the state space. The SSS met- 
ric is written in the form 



ds 2 



J-v 



dt 2 + N 2 dx 2 + r 2 (d0 2 + sin 2 6»d0 2 ) 



(1) 



where v, N and r are functions of the radial vari- 
able x. The function N is the radial gauge and below 
we will make a definite gauge choice which simplifies 
the field equations. The radial gauge function can be 
viewed as the analogue of the lapse function for spa- 
tially homogeneous models. Also, analogously, the 
expansion of the unit normals to the homogeneous 
hypersurfaces plays an important role in obtaining a 
compact state space. See M and references therein 
for a discussion of the lapse versus radial gauge and 
the use of the expansion for creating a compact state 
space. An expansion related function is given by 



6 = V3N~ 1 — (is + lnr) 
dx 



(2) 



variables according to 
dx 



K :=3r- 2 6" 



(3) 



A third dimensionless variable is defined as y := 
P/{p +P)- We shall also need the matter function 

F {y) '-=(p+p)-r = i-y- v^ 2 = i - v - 7 _1 , 

dp 

(4) 

where v 2 = dp/ dp, is the speed of sound and 7 := 
y~ 1 v 2 is the relativistic compressibility index. Now 
choosing N — V3yB _1 and using a prime to denote 
differentiation with respect to x, the field equations 
reduce to the dimensionless system 

£' = -yKT, + |P[1 + 2y(l - 2S)] , 
K' = 2y(l - K - 2P)K , (5) 

y' = -yF(y)E , 

where P := 1 — £ 2 — K, together with the decoupled 
equation for given by 



6' = -y(K - £ + 2£ 2 )9 



(6) 



The state space is defined by the inequalities |S| < 
1, P > and < y < yi < 1 where y\ is de- 
fined to be the lowest positive value of y satisfying 
F(y) = 0. The state space is a compact region hav- 
ing the form of a half pipe resembling a loaf of bread 
Q]. The bottom lies in the surface K = and the 
"roof ridge" stretches along the line {S = 0, K = 1}. 
The surface of the star is located on the curve de- 
fined by y = P = while the center corresponds to 
the points of the roof ridge 0. For the stiff fluid 
equation of state, y — p/(2p + p s ), implying that 
y varies in the range 0<y<2/i = l/2 where 



2/1 = Vo 



lira 



p — >oo 



Next we introduce two dimensionless gravitational 



y. Further F(y) = 1 - 2y so 
F tends to zero only in the limit of infinite pressure. 
The entire state space is therefore completely spec- 
ified in this case by the above range of y. Also, the 
system (o) is polynomial for this equation of state 
and hence regular throughout the state space. Note, 
however, that this nice and simple state of affairs 
does not always hold for more general equations of 
state for which y is not a monotone function of p. 
The function F(y) then has zeroes for values of y 
corresponding to finite pressures. This situation oc- 
curs e.g. for the Harrison- Wheeler equation of state. 
A more general formulation valid for a much larger 
class of equations of state will be given in . 



The compactness of an object has usually been 
measured by the single number R/M, i.e., the tenu- 
ity in the terminology of M. In the loaf of bread 
system the tenuity is simply given by 



R/M = 1 + 1/E. 



(7) 



However, as discussed above and in |jj, the tenuity 
by itself does not give the whole story for general 
equations of state. Instead the key to a more satis- 
factory understanding of stellar compactness is the 
effective potential which governs the null geodesies. 
It is also known as the centrifugal potential and is 
given by §] 



V = l{l + \)e lv r 



(8) 



A star is ultracompact precisely if V has an ex- 
ternum for some r > 0. The evolution equation 
for V takes the simple form 



V = -2y(l - T£)V 



(9) 



It follows that the extrema of V are located in the 
plane £ = 1/2. Note that while the right hand 
side of (f|) also vanishes for y — this is just a 
mathematical artifact caused by the vanishing of the 
chosen radial gauge at the surface. Of course, V may 
have an extremum at the surface but only if X = 1/2 
there. 

For any given equation of state there is a special 
"skeleton solution" representing the extreme case of 
infinite pressure at the center. Although unphysical, 
the skeleton solution is very important for under- 
standing the global properties of the solution space 
of SSS systems. In particular, the state space orbit 
of a typical relativistic configuration, after starting 
out from the roof ridge near y = y^ , spirals around 
the skeleton orbit before finishing at the surface, 
y = 0. It so happens that the skeleton solution 
for the stiff fluid can be given in the exact form 



£ = 



K = 



4(1 + y) 



(10) 



This is a special Tolman type V solution jl(j]. Thus 
the skeleton solution has V = and hence V — 
constant throughout its interior. It follows that 
there is a continuous family of circular null geodesies 
extending from the center to the surface of this stel- 
lar model. Now using (pi) and (JlTj) it follows that 
R/M = 3 for the skeleton solution. The model can 
be conveniently visualized by its optical geometry 
M which is shown in figure pi 





Figure 1: The equatorial submanifold of 
the optical geometry of the skeleton solu- 
tion for the stiff matter equation of state 
(p — p — p s ). The figure represents the star 
from r = 0.005 M to r = 9.5 M. The cylin- 
der is in reality infinitely long in the limit 
r — *• 0. The fat circles represent the surfaces 
r = 1, 2, . . . ,9 while the dotted curves repre- 
sent the surfaces r — 0.1 and r = 0.01. The 
shaded region represents the stellar interior. 
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Figure 2: Plot of £ versus Schwarzschild ra- 
dial coordinate for some selected values of 
the central density. The horizontal line at 
S = | represents the skeleton solution. The 
remaining curves correspond to the central 
densities (using clockwise enumeration about 
the upper left corner of the plot), p c /p s = 
1.329, 2.09, 6.2, 30, 274 and 9700. The most 
tenuous of the plotted models has R/M — 4 
while the next has R/M = 3 and is the first 
ultracompact model in the sequence. The 
third model has the smallest tenuity in the 
sequence, R/M = 2.75. It is a single neck 
model like the fourth and fifth models. The 
sixth model is a double neck model repre- 
sented by the dotted curve. Its inner neck is 
indicated by the circular disk while the posi- 
tions of its two bulges are given by the circles. 
Its outer neck is the Schwarzschild exterior 
neck at r — 3M. 



To understand the solution space it is necessary to 
analyze the full 3-dimensional system (0). However, 
for simplicity, in this communication we restrict at- 
tention to the behavior of E. The necks and bulges 
can then be read off as the zeros of the function 
£ — h. Necks and bulges correspond to unstable and 
stable circular null orbits respectively. The detailed 
geodesic structure requires a more careful analysis 



which is beyond the scope of this communication. 

The crucial observation comes from the realiza- 
tion that the orbits which correspond to sufficiently 
compact objects will spiral about the skeleton so- 
lution. This happens because the skeleton solution 
corresponds to a separatrix orbit going from a spiral 
critical point on the state space boundary at £ = h, 
K = y = i Since the skeleton solution sits on the 
surface £ = h the number of necks is determined by 
the number of revolutions in the Y,-K plane. In the 
limit of infinite central pressure the number of rev- 
olutions grows without limit and hence the number 
of necks can be arbitrarily large. In figure [2| £ is 
plotted as a function of the Schwarzschild variable 
r for some values of the central density including a 
value which leads to a double neck star. 

The null geodesic structure is intimately con- 
nected with quasi-normal perturbation modes of the 
relativistic stellar models || |j . These modes in turn 
are closely related to gravitational wave scattering. 
Threfore the appearance of multiple necks is likely 
to influence the picture of both gravitational pulsa- 
tion modes and gravitational wave scattering. The 
examples of multiple neck models we have given in 
this communication all correspond to unstable stel- 
lar objects. Although this would seem to limit the 
physical applicability of the results, moderately un- 
stable models may in fact have a role to play as 
intermediate stages in gravitational collapse situa- 
tions pjj . It is an open problem whether there ex- 
ist multiple neck models which are both stable and 
causal. We mention without proof that the stiff fluid 
equation of state admits a sequence of stable (and 
of course causal) single neck models but no stable 
double or higher neck models. 
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